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Abstract 

We strengthen a result from [T7j on the existence of effective interactions 
for discretised continuous-spin models. We also point out that such an 
interaction cannot exist at very low temperatures. Moreover, we compare 
two ways of discretising continuous-spin models, and show that, except for 
very low temperatures, they behave similarly in two dimensions. We also 
discuss some possibilities in higher dimensions. 
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1 Introduction 



If one tries to approximate a continuous-spin vector model such as the classical 
XY-model by a discrete approximation, whether for computational or for theoret- 
ical purposes [16] . one can in principle pursue two routes. 

1. Either one could consider the Gibbs measure for the original continuous-spin 
model, and discretise the spin by dividing the single-spin space into a large 
but finite number q of sets (intervals). By identifying all spins in such a 
set, one obtains a measure on a discrete-spin system. This measure can, 
but does not have to, be a Gibbs measure for some effective discrete-spin 
interaction. 

2. Or, alternatively, one could write down the same expression for the Hamilto- 
nian of the discrete-spin system as one has for the XY-model, and consider 
the appropriate Gibbs measure (s) for this discrete spin interaction. For the 
XY-model (the plane rotor) such models are called clock models, and their 
study goes back to Potts [21] . 

Here we present some results on discrete approximations of type 1, and com- 
pare them with what is known on discrete approximations of type 2. We will see 
that, except for very low temperatures, in two dimensions both approximations 
have quite similar properties, and we speculate on possible scenarios also in higher 
dimensions. 

We notice that a type 1 discretisation is appropriate for measuring purposes, 
and describes for example round-off errors, whereas a type 2 discretisation is what 
often happens in computer simulations. 

Generalisations to discretisations of more general continuous compact single- 
spin spaces are immediate, as follows from the analysis of [T7J [19]. 

Stated differently, and more formally, we can apply a local discretisation map 
T : S 1 i— >■ {1, . . . , q}, mapping a continuous local spin variable, taking values on 
the circle 8 , to its discretized image, before or after performing the Gibbsian 
modification with interaction $ relative to the product measure a. 

We then want to compare the images T\i of the Gibbs measures \i G G<f>, a 
of the initial model with a priori measure a and interaction $ with the Gibbs 
measures // G G<s>,Ta where Ta is the product of the a priori measures under 
the local coarse-graining T and G&,Ta are the Gibbs measures obtained from the 
specification which has the old interaction simply taken in the coarse-grained 
variables. 

We note, by the way, that conceptually such a question can be studied even 
more generally without making any assumptions on the first and/or second image- 
spin measure being discrete. 

A first important question to be asked is whether and when 

TG$,a = G$>,Ta 

for some 
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We remind the reader that if there is an $' such that T/i G Q^^Ta for a 
/i G £<j, iCn then TQ^^cQv ? a (3 • 

We are not aware that equality between the number of transformed Gibbs 
measures and the number of Gibbs measures for the transformed interaction - 
even if one exists- always holds, although we don't know of any counterexamples. 
We also note that, even under the assumption of equality of these two sets, one 
can have a different number of extremal Gibbs measures in the original and the 
transformed set, if different measures are mapped to the same one, as e.g. occurs 
in Fuzzy Potts models. 

We can also ask questions of closeness on the level of interactions, namely, 
what is the distance between original and transformed interaction <i($, $')? Fur- 
thermore, what is a good notion for the distance d here? 

There is the problem here that the spins, and hence the interactions, live on 
different spaces, one discrete, one continuous. If one compares the two discretisa- 
tions one has at least the advantage that the corresponding interactions will live 
on the same space. 

What can be said about closeness of the measures fi G to jl G ^^q? This 
question is subtle, since we could look here for closeness on local observables,-that 
is, in the weak topology-, on the level of long-range characteristics like decay of 
correlations, on the level of the phase diagram in parameter space, ... Here we will 
call the two discretisations close, if the two models have a similar phase diagram 
and/ or similar correlation decay in their Gibbs measures. 

We note that the proof concerning the locality properties in [17] for type-1 
discretisations makes essential use of the Dobrushin uniqueness theorem, even 
though we need not be in the uniqueness regime and in fact are allowed to be in 
a phase-transition region when we discretise. This will also be the case here. 

Discretisations can be viewed as single-site coarse-grainings, similarly to the 
fuzzification or amalgamation of discrete-spin systems as treated in e.g. [21 [131 EH 
[28] . but now the "fuzzification" goes from a continuous "alphabet" to a discrete 
one. 

2 Gibbsianness of discrete approximations of the 
XY lattice model 

2.1 Notation and Definitions 

We will consider lattice spin systems with a single-spin space Sl , on a lattice Z d , 
and a configuration space Vt = Q Z . We will mainly consider the XY-model, 
for which Q Q is the circle S 1 , and discrete approximations thereof, in which S 1 is 
divided into q equal arcs of length We will indicate the spin variables at site i 
(which always will be elements of the unit circle) by <7j, Ui, and similarly spin 
configurations in a volume A by oa, <^a, Va- 

We will consider Gibbs measures, which are defined for (here translation- 
invariant) absolutely summable interactions $ (that is, YIa- oga I I < °°) v ^ a 
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the DLR equations, expressing that given an external configuration rj A c, the prob- 
ability density of configurations in a volume A is given by the Gibbs expression 



da A 



exp (-H A (a A i] A c)) 



Z A 



where H A (a A rj\c 



2j (3^ A (a A r] A c), 
A; AnA^0 

(1) 

and a A is the product of a over the sites in A. This should hold for all volumes 
A, internal configurations a A and external configurations r] A c The corresponding 
collection of (everywhere instead of almost everywhere with respect to the Gibbs 
measure defined) conditional probabilities forms a "specification" . See e.g. [U [7J 
T2] . In the standard nearest-neighbour models, (the plane rotor or XY-model), as 
well as in the clock models, where the spins take discrete values, we have 



-H A (a A r] A c) = (3 ^2 



Oi ■ a A 



p E 



(2) 



<i,j>€A 



<ieAjeA c > 



2.2 Conservation of Gibbsianness under local transforma- 
tions: fine discretisations 

One of the main results of [IT] , see also [19] , concerns conditions under which a 
discretisation of a Gibbs measure is again Gibbsian. These results were obtained 
as corollaries to a theorem on the preservation of Gibbsianness which also holds 
for much more general types of local transforms like time evolutions. So, it is 
worthwhile to reconsider specifically the local transformations. 

The condition for preservation of Gibbsianness is temperature-dependent, and 
the main example we want to discuss here is the discretisation of the circle into q 
equal arcs. At inverse temperature /3 the result implies that for q, dependent on 
0, large enough, the discretised measure is a Gibbs measure. 

To be more precise, suppose for each I 6 S' :— {1,2,- •• ,q} we denote by 
S>j the Zth arc of the circle S 1 cut out by the discretisation operator T. Then, 
Tail) = ai§j). Given // e G<f> >a , one of the main results in [17], Theorem 2.5, is 
that the discretised measure T\i is Gibbs if 



sup 



E 

3& d \{i} 



Cij < 1, 



(3) 



where 



2C 



■1 j 



*up„_.. r . . / s i a(d<Ji 



'3' '3 

ies 



if \i 



otherwise. 



(4) 



Thus, constrained on a discrete-spin configuration, the constrained system must 
be in the Dobrushin uniqueness regime uniformly in the chosen constraint. 
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While looking for good upper bounds for the right-hand side, we can at not 
much additional cost revisit the more general situation and give an improvement 
to the criterion from [T7] for Gibbsianness for local discretisations. 

We put ourselves in a slightly more general context than that of the discretisa- 
tions in [T7], and we will take the local spin space S just to be a general compact 
measurable space. No a priori metric is given; it will be produced by the Hamil- 
tonian itself. As in [17] , let a decomposition be given of the form S = Us'es' ^ s '- 
Here S' may be a finite or infinite set. Put T(s) := s' for S s > 3 s. This defines a 
deterministic transformation on S, called the fuzzy map (the discretisation). 

Now we deviate from [17]. Let G be the vertex set of a general graph and 
define a family of metrics (dy)jeG\{i} on t ne local spin space at the site i G G by 



where for any i e G, i c = G\ {i}. It is important here (as well as in the formula 
specific to the rotors above) that the supremum is taken over spins (, ( which are 
constrained to take the same coarse-grained image at j. We are allowed to do 
this since we are analyzing the constrained system. In this way the metric at the 
site i depends also on the size of the coarse-graining at j. The metric measures 
how strongly a variation at the site j can maximally change the difference in 
interaction energy between local spins cr^Tj. 

Our criterion of the fineness of the decomposition will involve the correspond- 
ing j -diameter, namely the quantity diamy(A) = sup s teA dij (s, t) where A runs 
over the sets in the decomposition. 

Theorem 2.1 Let \i be a Gibbs measure of the specification with Gibbsian poten- 
tial $ with an arbitrary a priori measure a, on a graph with vertex set G. Let 
T denote the local coarse- graining map where we assume that a(S s >) > for all 
labels s' £ S". 
Suppose that 



Then the transformed measure T(fi) is Gibbs for a specification 7' with an 
absolutely summable discrete-spin interaction 

In all cases this is an improvement over the criterion of [T7j (which we don't 
repeat here in detail, because it requires the introduction of additional structure 
which we don't need here.) 

It is also an improvement over what a direct application of the high-temperature 
version found in Georgii [12] would give for our constrained model. That would 
only give a bound in terms of the right-hand side of the inequality of the form 




sup 

C jC =C jC ;T(C i )=T(c i ) 





(6) 




(7) 
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Here we discuss the application to rotor models. Consider first the rotor model 
on a circle S 1 . We have for n.n. i and j by Cauchy-Schwartz that 

dijipu n) = f3 sup {(Ti - n) ■ (Q - Q) 

< o-j - n 2 2sin- 
Q 

and so dianijjS' s / = (3 x (2sin^) 2 . This gives the criterion 

2c^(sin-) 2 <l (9) 
Q 

for Gibbsianness of the coarse-grained model. Note that the standard estimate 
(J7J would give a worse condition without the square. 

For a local spin space which is a d-dimensional sphere not much changes. The 
formula for the metric dij stays the same. Let us assume that ip is one half of 
the maximal angle under which a set S s > appears as seen from the origin. This 
quantity is a measure of fineness of the discretisation. Then, going through the 
same steps, we obtain as a criterion for Gibbsianness that 

2d/3(sin^) 2 < 1. (10) 

Proof of the Theorem: The proof follows as in [T7], by estimating CV,. This 
constant is a bound on the Dobrushin interaction matrix of the initial model con- 
ditional on the transformed spins, uniformly in the values of the transformed spins. 
In particular, for each site ieG, CV, is a uniform upper bound on the variational 
distance between the "first-layer models" in {i}, in which <ji is constrained to take 
values in S s > for some fixed prescription of partitions given by the image spins a' { , 
wrt external configurations that coincide everywhere except at site j. 

More precisely, we take two conditioning configurations in the original (first- 
layer) model, (,rj G Q with Qc = r]jc and denote by u (ai) = —H^aiQc) and 
M i( cr t) = —Hi{ciVi c ) the corresponding values of the single-site Hamiltonians an- 
chored at i. Defining u t = tu\ + (1 — t)u , h t = e Ut ls ,/a(e Ut ls ,) and \ t (dai) = 

h t (ai)a(d<Ji), with t e [0,1], we note that \ (d<Ji) = ^^da^Q and Xi(dai) = 
■ji(dai\r]), where 7j's are the single-site parts of the conditional distributions (ker- 
nels) of the initial model obtained via ([T]) after replacing 5" with S a > . Note however 
that the constraining configuration o' does not appear any more in the notation, 
for the sake of simplicity. 

Now comes the estimate which for coarse-grainings improves the one from [17] 
(in which, however, also general transformations beyond coarse-grainings were 
treated), and also ([7J from Georgii. For the first step of the proof, we obtain the 
following bound for the total variational norm of Ao — Ai; 

1 , d , , 
dt—h t {ai 
dt 



2||A Q - Ai|| = J a(doi)\hi(ai) - h (<Ti)\ = J a(da { 
< J dt J Xt(dai) Hi(aiCic) - ^(cr^c) - J \ t (dTi) (H^nQc) - B. i {r i r] i . 
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The simple but essential next estimate will be uniform in the image measure of 
X t , under >->• H^aiQc) — Hifarjic) =: /^(crj). Namely, by further making use 
of the notion of the j-diameter of the set S a / to bound the variation of the energy 
terms we get 



f^(X t )(dx) / f<«(X t )(dy)\x-y\ 



(12) 



< sup / X(dx) 



dyX(dy)\x-y\ 



D 



where D = diam^ S a ' / 2 and the supremum is over the probability measures A 
on the interval [— D,D\. For this supremum we use the following lemma (after 
scaling with D). 

Lemma 2.2 For all probability measures p on [—1, 1] we have 
Q(p) := / p(dx) f p(dy)\x — y\ < 1 with equality for p = |(5i + 

Note the improvement over the simple upper bound 2. Observe also that 
the upper bound on 2||Ao — Ai|| (TTTT) obtained via f|T2l is independent of rj, C G 
Q and constraint a\. Therefore f|T2|) provides a uniform upper bound on 2CV,-. 
By scaling up the interval in the Lemma with a factor D and putting together 
our previous estimates we get 2CV,- < sup s / diam^ S a >./2, and hence Ysj-.j&C'ij — 
\ ^2jeG\i su Ps' diam i j(S'< j /). The rest of the proof follows from the definition of the 
Dobrushin constant. This proves the Theorem. □ 



For the sake of completeness let us also give an elementary proof of the Lemma. 

Proof of the Lemma. By density arguments we can approximate any p by 
convex combinations of finitely many Dirac measures of the form Y^=i Pi^x t where 

X{ ^ Xi-\-\. 

Let us look at Q as a function of the Z-th location, keeping the other locations 
fixed, and keeping the p^'s fixed, xi h> Q(Y^=xPi^xi)i where X[ is constrained 
to be greater than or equal to its left neighbor and less than or equal to 
its right neighbor This function is linear. Hence the function takes its 

maximum when x\ becomes equal to one of its neighbors. This shows that the 
maximum of Q over the set of combinations of n Dirac measures is dominated 
by that over combinations of n — 1 Dirac measures. Iterating this argument we 
see that the maximum of Q over all probability measures is reached for a linear 
combination of two Dirac measures p5 x + (1 — p)5 y . Noting finally that the max 
over Q(p5 x + (1 — p)5 y ) = 2p(l — p)\x — y\ is reached for p we are done. □ 
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As far as the bound on X^ez d \{i} A? * s concerned, the above result is an 
improvement over Theorem 2.9 of [17] which however was formulated in a much 
more general situation. Indeed, the latter gave rise to the bound Yujez d \{i] — 
4d *P eP . The more general set-up of [T7J allowed also to treat (partially) stochastic 
single-site maps, such as infinite-temperature stochastic dynamics. The estimates 
on the Dobrushin constant used there were of the "high-field" type , whereas here 
we make use of a "high-temperature" version. 

Once the refinement is large enough (q very large at a fixed temperature) the 
effective interaction has as its dominant term the nearest-neighbor interaction of 
the clock model. 

We notice that the discretised model inherits various properties from the orig- 
inal XY-model. In particular, if the correlation functions decay slowly, as they do 
in two dimensions at low temperature when one is in a Kosterlitz-Thouless phase 
this remains true after the discretisation. 

The continuous symmetry of the original model is also inherited. In three and 
more dimensions there is a continuum of Gibbs measures, as the circle symmetry 
of the original XY-model is broken. Discretising the spin space results again in 
a continuum of Gibbs measures (which now however are not related to a broken 
symmetry of the discrete-spin model). 

We have thus proved the following theorem. 

Theorem 2.3 For each d > 3 there is a such that for q > q there is an 
interaction $' with a discrete -clock - rotation invariance such that there are un- 
countably many translation-invariant ergodic states in the set of Gibbs measures 
Q<$>i (taken with uniform a-priori measures). 

This argument provides an independent rigorous route to the existence of 
an intermediate enhanced-symmetry Kosterlitz-Thouless phase in a discrete-spin 
model, combining our general criteria for preservation of Gibbsianness under local 
coarse-grainings with properties of the original continuous-spin model. 

3 Comparing the discretisations 

At high temperatures, in the paramagnetic regime, everything is well-behaved, 
but not of great physical interest. We will therefore discuss what happens in 
subcritical-temperature regimes. 

It is a remarkable fact that the standard nearest-neighbor large-g clock model 
in two dimensions has the property that there is a Kosterlitz-Thouless phase 
with slow decay and an enhanced continuous symmetry occurs at an intermediate 
temperature regime [SJE]. On symmetry enhancement, see also [T5] . 

The values for which this occurs are such that q should be large enough, for 
a given low temperature. As we have just seen, the discretised XY-model can be 
described by a summable interaction -in which the nearest-neighbor terms are the 
dominant ones-, in just such an intermediate regime. 
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On the other hand, the nearest-neighbor clock model at fixed q and at very 
low temperatures (/3 > 0(q 2 )) will have q ordered phases, that is q different Gibbs 
measures, similar to the q ground states, all with exponential correlation decay. 
This follows directly from a Pirogov-Sinai argument. On the contrary, we can 
show that for a fixed even q once the temperature is low enough (how low depends 
on q), the -type 1- discretised Gibbs measure becomes non-Gibbsian. Indeed, if 
we take an alternating configuration for the discretised spin this implies that 
alternatingly the spin is either in the most Northern (on sites in one sublattice) 
or the most Southern interval (on the other sublattice) of size 2zl. We argue 
that such a configuration is a point of essential discontinuity for a conditional 
probability of the discretised measure. Conditioned on this, the original spins 
(which are forced by the constraint on which we condition to be almost opposite, 
but by their interaction prefer to be pointing in the same direction), will have two 
ground states, one pointing alternatingly North- West, South- West, and the other 
one alternatingly North-East and South-East. The deviations in the Western , cq 
Eastern, direction are of order O(^), which means that the energy gap between the 
two ground states is of order O(^). Therefore, at sufficiently low temperatures 
(j3 > 0(q 2 )), there will be two different Gibbs measures for the constrained model, 
and this will imply the non-Gibbsianness of the discretised measure. The details 
of the argument can be worked out in a straightforward manner along the lines of 
El 122], see also [3]. 

Thus the analogy between the two discretisations breaks down just in this 
very-low-temperature regime. The measures then are not even close any more 
on the level of local observables. Since in this regime one finds very different 
behaviour, one discretisation resulting in a non-Gibbsian measure, and the other 
one in q different Gibbs measures. 

In higher dimensions, for the XY-model there is a continuum of Gibbs measures 
[TT] at low temperatures, which, as indicated above, are mapped to a continuum 
of different Gibbs measures for the discrete spins in an intermediate regime. 

It would be interesting to see if the restoration of continuous symmetries which 
happens for the two-dimensional clock model would have a higher-dimensional 
analogue, in that in some intermediate-temperature regime there might exist a 
continuum of Gibbs measures, even for the nearest-neighbor clock model. We 
conjecture that the intermediate phase studied in [T5J [201 1221 I2U I2E1 I2Z], might 
be of this type. In the terminology of Ueno et al [2Z] we would have a continuum 
of "Incompletely Ordered Phases" , where the order can be in the two spin direc- 
tions n, n + 1, where n mod q e {1, • • • , q}, with continuously varying weights of 
these directions. Although there seems to be some doubt whether there exists an 
intermediate phase at all in a region in between the q ordered, ground-state-like, 
phases and the high-temperature paramagnetic phase, the numerical results up 
to now for the nearest-neighbor clock model appear to be inconclusive. It there- 
fore seems worthwhile to investigate if an "enhanced-broken-symmetry phase" as 
decribed above, which can be obtained by discretizing a continuous- spin model, 
could also occur for the nearest-neighbor clock model. 

The breakdown of the analogy at very low temperatures holds for the same 
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reason as in two dimensions. The only property we used was the bipartiteness 
of the lattice. If we choose the North in the direction of the magnetisation, the 
arguments are unchanged. 

As we have just seen the transition between the Gibbsian behaviour and the 
non-Gibbsian behaviour occurs at > 0(q 2 ) for type-1 discretisations. 

The analysis of [8] for type-2 discretisations similarly appears to provide a 
transition value at /3 — 0(q 2 ) for the transition. For some numerical results, 
indicating this asymptotics in more detail, see [25] . 

A heuristic reason for this behavior is that the model in the scaled variables 
q times spin-angles with discretisation width 1 approximates a discrete Gaussian 
model at effective inverse temperature (3/q 2 ; if this parameter is below the value for 
the roughening transition (which is rigorously known to take place in the discrete 
Gaussian) the model behaves like a massless Gaussian, while above it behaves like 
a massive model in the Peierls regime. Compare also Theorem C on page 40 of 
[TO] . There it is mentioned as a conjecture that, at fixed q > 5, the threshold 
values in temperature between low temperature regime and intermediate regime 
on the one hand, and intermediate and high temperature regime on the other 
hand, should be sharp, and different. See also [1] for some numerical support for 
this. 

4 Conclusions 

We showed how to compare two different ways of discretising spin models, namely 
either starting from the Gibbs measures, for which we have a controlled approx- 
imation in a temperature-dependent regime, or starting from the interactions. 
We extended the regime in which we have such a controlled approximation of 
the discretised Gibbs measure, and also pointed out that it cannot be extended 
to very low temperatures. Thus the results are essentially optimal. In the two- 
dimensional XY-model both discretisations display the same Kosterlitz-Thouless 
phase in an intermediate-temperature regime. 

As for the higher-dimensional case, we suggested the possibility of an enhanced- 
continuous-symmetry-breaking phase occurring at a region of intermediate tem- 
peratures in discrete <i-dimensional clock models for d at least three. 
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